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Fine structure in the gamma-ray sky and the origin of UHECR 

A.D.Erlykin a *, A.W.Wolfendale b 
a P.N.Lebedev Physical Institute, Moscow, Russia 
b University of Durham, Durham, UK 

The EGRET results for gamma ray intensities in and near the Galactic Plane have been analysed in some detail. 
Attention has been concentrated on energies above 1 GeV and the individual intensities in a 4° longitude bin 
have been determined and compared with the large scale mean found from a nine-degree polynomial fit. 

Comparison has been made of the observed standard deviation for the ratio of these intensities with that 
expected from variants of our model. The basic model adopts cosmic ray origin from supernova remnants, 
the particles then diffusing through the Galaxy with our usual 'anomalous diffusion'. The variants involve the 
clustering of SN, a frequency distribution for supernova explosion energies, and 'normal', rather than 'anomalous' 
diffusion. 

It is found that for supernovae of unique energy, and our usual anomalous diffusion, clustering is necessary, 
particularly in the Inner Galaxy. An alternative, and preferred, situation is to adopt the model with a frequency 
distribution of supernova energies. The results for the Outer Galaxy are such that no clustering is required. 

If their explosion energies are distributed then supernovae can be the origin of UHECR. 



1. INTRODUCTION 

Although supernova remnants (SNR) are often 
invoked as the source of cosmic rays (CR), of en- 
ergy below the 'knee' at about 3 PeV, [l|2j , there 
are still a number of imponderables. Here, we 
concentrate on two aspects: 

(i) SN energies are not unique, but have a cer- 
tain frequency distribution involving some 'super- 
supernovae' [3]. We refer to this model as SSNR. 

(ii) SNR in the Galaxy are not distributed at ran- 
dom, but rather in clusters, in both space and 
time and, preferentially, in spiral arms [4]. This 
is NSNR model ( Normal SNR ). 

In our previous work [5] we used the accelera- 
tion model of [2J, which assumed that particles 
are trapped for up to 8.10 4 y after the SN ex- 
plosion and that SN occurred randomly in time 
in the Galaxy (at an average rate of 1CU 2 y ). 
The distribution in space was drawn from the 
usual average radially symmetrical distribution 
of SN surface density, with a peak at a Galac- 
tocentric distance of D=4 kpc and falling rapidly 
with increasing D in the vicinity of the sun ( at 
D=8.5 kpc ). Here, we go on to make allowance 
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for the clustering of the SN - and hence the SNR 
- in various ways, and, separately, consider a fre- 
quency distribution of SN energies: the SSNR 
model. 

The quantity used as an indicator of CR inten- 
sity in the Galaxy is the gamma ray intensity [5] 
which has been measured as a function of longi- 
tude, latitude and gamma ray energy. Most of 
the present work relates to gamma rays above 1 
GeV which are considered to be produced pre- 
dominantly by CR nuclei (mainly protons) of me- 
dian energy ~ 40 GeV [7]. In zero'th order of 
accuracy, the mean CR intensity along a partic- 
ular line of sight (£, b) is simply proportional to 
the measured gamma ray intensity divided by the 
column density of target gas. In the actual analy- 
sis we use a more elaborate method but the basis 
is the same. 

The method is to study the ratio of the ob- 
served gamma ray intensity to the large scale av- 
erage and use that as an indicator of the small 
scale variability of CR intensity from place to 
place. The variability can then be compared with 
that predicted by models involving SN, clustered 
and unclustered, normal and with different ener- 
gies, using the same technique of comparing the 



1 



2 



intensity with the average. 

2. THE BASIC DATA 

Figure 1 gives an example of the data ( from 
[6] ); it refers to gamma rays with energy above 
1 GeV and having \b\ < 2°, i.e., it relates to the 
Galactic Plane. Identified discrete sources have 
been removed and the workers estimate that less 
than 10% of the remaining flux is due to unre- 
solved sources. 
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Fi gure 1 . Profile of the gamma ray intensity vs longi- 
tude, for E 1 > 1 GeV and |6| < 2°. The results are from 
[6]. The standard deviations on the points due to Poisso- 
nian fluctuations are typically ±2.3% in the Inner Galaxy 
and ±4.8% in the Outer Galaxy. The experimental data 
are averaged over the longitude bin Al = 1°. The smooth 
line is the best fit of the profile by the 9-degree polynomial. 



We have used the basic data for the numbers 
of gamma rays to evaluate the statistical uncer- 
tainty in the individual intensities. 

The data have been divided into 'Inner' {£ : 
270° - 350°; 10° - 90°) and 'Outer' (£ : 90° - 
180° - 270°), the innermost region (\£\ < 10°) 
being omitted because of confusion there. 



Two latitude ranges are considered: \b\ < 2.5° 
and \b\ : 2.5°/5°. The gamma rays from the indi- 
vidual b- and /-ranges are derived from regions of 
the Galaxy distributed along the line of sight but 
we can identify a median D- value for each range. 
Restricting attention to gamma ray production 
in molecular clouds with the scale height of ~ 70 
pc, the corresponding distances are ~3 kpc and 1 
kpc from the sun, for \b\ < 2.5° and \b\ = 2.5°/5° 
respectively. Allowing for production in HI, with 
its greater scale height, the D-values will be some- 
what different ( but by not much because more 
of the gas - associated fluctuations arise from the 
clumpy Hi rather than HI ). 

3. THE ANALYSIS 

3.1. The Measured Fluctuations 

The 'figure of merit' adopted is the ratio of the 
observed gamma ray intensity to the smoothed 
line through the points ( 9-degree polynomial, 
equivalent, approximately, to a smoothing over 
a 20° longitudinal bin for the Inner and Outer 
Galaxy regions, separately ). The shape of the 
line is determined essentially by the kiloparsec- 
scale spatial variations of the cosmic ray intensity 
and the target gas. Figure 1 gives the smoothed 
curve for |6| < 2°. 

The ratios have been determined for various 
bins of longitude, M, specifically, 1°, 2°, 4° 
and 8°. Figure 2 shows the corresponding fre- 
quency distributions of the differences between 
the EGRET data and the fit. The numbers of 
entries, the means and the rms (cr— ) values are 
given in Table 1. It is with the latter rms values 
that we are mainly concerned. 

The appropriate tr-value to be used is a com- 
promise between the A£ over which intensity 
changes are expected (all values but with a bias 
towards degree scales) and the need to get away 
from small A£-values, where bin-to-bin correla- 
tions are serious, because of the finite range of 
latitudes (\b\ < 2.5°) over which the averaging is 
made. 

We consider that A£ = 4° is most appropriate 
and this will be used in the model calculations. 

Figure 3 gives the values of a (4°) for the ob- 
served data ( after correction for noise ) and 
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Figure 2. The distribution of differences between the 
EGRET intensities and the polynomial fit for various /\£- 
values. The numbers of entries, the means and the RMS 
(it—) values are given in Table 1. 



for various model calculations, to be described 
shortly. The values are plotted at increasing 
(from left to right) distances from the Galactic 
Centre of the median gamma ray production re- 
gions. 



3.2. Variations in the column density of 
gas 

In all our models the radial distribution of gas 
has axial symmetry. It means that for a smooth 
non-fluctuating CR distribution the longitudinal 
distribution of the gamma-ray intensity has also 
to be smooth. However, we know that the ISM is 
highly non-uniform and this non-uniformity will 
contribute to the fine structure of the gamma-ray 
intensity profile. 

A number of facts are relevant, as follows: 

(i) The large scale dependence of N(H) on I is 
automatically included through the radial distri- 
bution. 

(ii) The 9-degree polynomial fit causes changes 
on scales of about 61 — 20° ( and above ) to be 
allowed for. 

(iii) Much of the dumpiness of the gas in the ISM 
comes from molecular clouds and we allow for cor- 
relation between the SNR and the clouds in the 
analysis. 

(iv) Results of the type shown in Figure 2 show a 
dependence of a on Al in the range Al : 1° — 8° 
which is very similar for all the models and the 
observed values. 

Although the implication of the above is that 
the effect of the fluctuations in N(H) ( both 
atomic and molecular ) is not likely to be dra- 
matic, it is certainly finite. Its magnitude will 
be greatest in the Inner Galaxy at low latitudes, 
\b\ < 2.5°, where contributions from the 'far' In- 
ner Galaxy are important. 

Estimate for both the Inner and Outer regions 
have been made. For N(HI) the estimated value 
of the logarithmic standard deviation is 0.027 us- 
ing the well known column densities of HI [5]. 
For N(H 2 ), data [9] give a value 0.026 referred 
to the whole column density ( N(HI) + N{H2) ), 



leading to an overall value of a g 



0.037. Com- 



parison can be made with the observed value of 
a bs for this l,b region: 0.080. Subtraction in 
quadrature leads to a = 0.071, i.e. a small, but 
finite reduction. 

Support for our analysis comes from a study 
of the EGRET observations and analysis pub- 
lished prior to the latest observations used by 
us [816] . In this work, predictions were made 
for a model in which the column densities of gas 
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were used together with an inferred CR inten- 
sity distribution which was correlated with the 
volume density of gas ( derived from the col- 
umn density and assumed rotation curve ) but 
'smoothed' spatially at the kpc level. The major- 
ity of the excursions in predicted intensity were 
thus due to the gas column density variations 
alone. The fluctuations of the observed intensi- 
ties ( E 1 > lGeV, \b\ < 2°, \l\ : ±90° ) about the 
prediction have a = 0.06 ± 0.015. This value is 
quite consistent with ours. 

Turning to regions away from the Galactic 
Plane in the Inner Galaxy ( |b| = 2.5° - 5° ) the 
corresponding value for HI and H2 is cr gas = 0.04. 
In the Outer Galaxy the values are a = 0.025 for 
|6| < 2.5° and 0.035 for |6| = 2.5°/5°. It will be 
noticed that in all cases the corrections will be 
small. 

Our contention is that the major cause of 
the fluctuations in gamma-ray intensity on the 
4° scale is the stochastic nature of the SNR in 
space and time and, probably, in energy output. 
The column density of gas undoubtedly fluctu- 
ates with I but the observed gamma-ray intensity 
fluctuations are much greater because of the CR 
intensity variations along the line of sight. 

3.3. The model predictions 

Our 'usual' model |2|5j comprises CR produc- 
tion by SNR, as described in paragraph 1, and 
propagation with 'anomalous diffusion'. The 
SNR were assumed to be all identical and inde- 
pendent in time and space ( i.e. NSNR ), within 
the constraints of the dependence of the mean 
SNR density on Galactocentric radius. We de- 
note this as NC ( non-clustering ) , and the results 
for this situation are shown in Figure 3. This is 
the NSNR model, but with no clustering. 

The SSNR model has the same mean energy 
content per SNR as in the previous case but a 
distribution of energies about the mean, with a 
frequence distribution of log energy having a log 
standard deviation of 1.2, following [3 . The re- 
sults for this type of super-supernovae are now 
indicated as SSNR in Figure 3. 

Turning to the possibility of the clustering of 
NSN, our clustering model has the following ide- 
alised form. At any one place there are SN in sets 
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Figure 3. Summary of standard deviations (cr- values) of 
gamma ray intensities about the best-fit: (a) - normal dif- 
fusion, (b) - anomalous diffusion. Small downward correc- 
tions have been applied to the OBS ( observed )-values to 
allow for experimental noise. Model simulations are: NC 
- non-clustering, C - clustering, SSNR - with a frequency 
distribution of energies. Mean standard deviations are ob- 
tained by averaging 4-8 independent samples. Errors are 
not indicated for clarity of the figures, but the distribution 
of cr for different samples is broad due to the stochastic 
nature of SN explosions. Typical errors for anomalous dif- 
fusion are 0.002-(NC), 0.0015-(C) and 0.014-(SSNR) for 
the Inner Galaxy, 0.004-(NC), 0.004-(C) and 0.02-(SSNR) 
for the Outer Galaxy. For normal diffusion the errors are 
by 2-5 times less. The Galactocentric distances are the 
means from which much of the gamma ray intensities are 
derived. 



of 10, distributed at random, in time, over 10 6 y. 
The results for this model 'C ( clustered ), are 
also shown in Figure 3. Here, the SN are identi- 
cal. 

Finally, calculations have been made for the 
same situation, as indicated, for 'normal diffu- 
sion' (ND) rather than our standard: 'anomalous 
diffusion'. 
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3.4. Comparison of observations with the 
models 

3.4.1. General remarks 

It is apparent that the observed a values are 
higher than expected for the NSN model without 
clustering - i.e. non-clustering (NC). For normal 
diffusion NC(ND) the discrepancy is bigger still. 
Below we examine various factors which influence 
our model predictions. 

3.4.2. Mixed primary mass composition 

The model predictions shown in Figure 3 re- 
late to primary protons. However, the observed 
gamma rays are produced by all CR nuclei, i.e. 
by the mixed primary mass composition. Primary 
nuclei are more efficient in producing gamma rays 
than protons of the same energy per particle due 
to the higher interaction cross-section and higher 
number of interacting nuclcons in the collision. 
At the same time fluctuations of the gamma-ray 
intensity created by CR nuclei are lower than 
those for protons, therefore our model predictions 
should be reduced for the mixed primary mass 
composition. We calculated the fluctuations for 
the mixed composition with relative abundances 
equal to 0.416 for protons, 0.265 for helium, 0.170 
for CNO group, 0.149 for iron group of nuclei 
[TP] , The result is that we have to multiply our 
predicted cr-values by 0.93, 0.86, 0.85, 0.88 re- 
spectively for the increasing galactocentric dis- 
tances indicated in Figure 3. The uncertainty 
of this correction is within the range of 0.005 - 
0.019. As a result of the correction the discrep- 
ancy between observed and expected cr-values for 
the non-clustering case (NC) with a mixed pri- 
mary mass composition is bigger still. 

3.4.3. Mode of diffusion 

It is relevant to examine the question of 'normal 
or anomalous diffusion ?'. The mode of diffusion 
has obvious implications for the present work. 
Calculations have been made for all cases, i.e. 
non-clustering, clustering and SSNR models us- 
ing normal diffusion, instead of anomalous. The 
result is that the fluctuations and cr-values are 
lower for normal diffusion in all cases ( see Figure 
3a compared with Figure 3b ) , thus the difference 
between the non-clustering case and observations 



will be even larger. 

3.4.4. Super-supernovae 

Examination of Figure 3 shows that the cr- 
values for SSNR with a distribution of explosion 
energies is higher than for clustering (C) and for 
the observed fluctuations (OBS). It is due to the 
basic assumption about the fluctuations of the ex- 
plosion energy. Added to the fluctuations caused 
by the stochastic distribution of SN in space and 
time these fluctuations result in a rise of the total 
fluctuations. Inevitably, it is possible to 'dilute' 
the SSNR values by assuming that only a fraction 
of the SNR are of such high energy. According to 
[3] the fraction of such SSNR is equal to 0.2 for 
SNIbc and 0.1 for SNIIn, but we regard it as a 
quantity 'to be determined'. 

3.4.5. Clustering of supernovae and the 
presence of SSNR 

Although at first sight the case for some clus- 
tering of SNR at all distances is strong there is 
a complication that must be allowed for. This 
concerns the correlation of SNR with target gas, 
particularly in molecular form ( principally H2 ) . 
This correlation is on linear scales smaller than 
the kpc-scale already allowed for. Our 4° bin of 
longitude would correspond to a linear scale of 
280 pc at a distance of 4 kpc ( the typical distance 
to the production region in the Inner Galaxy for 
|b| < 2.5° ). Such a distance would embrace likely 
distances of molecular clouds (MC) from some, at 
least, of the SNR. 

The magnitude of this correlation is not clear 
but some progress can be made. The list of known 
SNR [11] has only 3% of the entries mention- 
ing adjacent MC but this fraction will increase 
with increasing SNR-MC distance. Inspection of 
the Galactic maps of MC [9] in the Inner Galaxy 
shows a mean separation of about 500pc; thus a 
correlation factor of about 0.5 would appear to be 
indicated. For the general case, it appears that 
the relation between the various quantities for the 
situation where there is partial clustering of the 
SN ( coefficient / ) and a partial correlation be- 
tween SNR and molecular gas ( coefficient 0.5 ) 
is: 
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where a b s is the total, observed (corrected) stan- 
dard deviation, a gas characterizes the fluctua- 
tions of the column density of the gas ( as deived 
in §3.2 ), a nc , a c and a ssnr are the expectations 
for non-clustering, clustering and SSNR models 
respectively and the mixed primary mass compo- 
sition, R nc , Rc and R ssn r - ratios of the model to 
the observed gamma-ray intensities for the same 
three models, g is the fraction of gas in molec- 
ular form and 5 is the fraction of SSNR among 
the total SNR. The first term on the rhs of equa- 
tion (1) relates to SNR in the ordinary inter- 
stellar medium, the second term - to SNR in- 
side MC. The equation satisfies the lower limit 
<?obs = <?nc when / = S = g = 

In the calculations it is assumed that SSNR 
are so rare that they are not clustered, whereas 
half of the ordinary SNR (NC) are clustered in 
MC. The values of g are, for the approximate me- 
dian D-values relating to Figure 3: 0.5, 0.3, 0.1 
and 0.1. In fact, different models give different 
absolute gamma-ray intensities, but they can be 
easily normalised to the experimental data mul- 
tiplying the explosion energy by a constant. The 
value of relative fluctuations a in this operation 
remains unchanged. Therefore, we adopted all 

nc Rc — Rssnr 1- 

The fraction / is therefore the only unknown 
variable and it can be calculated from equation 
(1). The results are shown in Figure 4. 

4. DISCUSSION 

The situations for the two models can be con- 
sidered in turn. 

For the NSNR model, the model used by us 
in much of our work, considerable clustering is 
needed in the Inner Galaxy. Indeed, with / > 1 
the clustering must exceed that adopted by us, 
viz. 10 SN per I0 6 y coincident in position. Al- 
though not impossible, it is unlikely that the nec- 
essary higher number ( 20-30 ? ) is present, even 
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Figure 4. The fraction of the clustered SNR as a function 
of galactoccntric distance for normal ( o ) and anomalous 
( * ) diffusion and for 3 SSNR fractions: (a) - <5 = 0, (b) 
- <5 = 0.15, (c) - a = 0.2. The lines connecting the points 
are 3-degree polynomial fits to the data. The abscissae 
for the two curves are slightly displaced with respect to 
each other to avoid an overlap of the errors. For 5 = 0, 
the first point ( extreme left ) for anomalous diffusion is 
at / = 3.21. 



in the congested region of the Giant Molecular 
Ring at D « 4fcpc. 

Instead, we incline towards the SSNR model. 
Inspection of Figures 4b and 4c shows that in 
the important Inner Galaxy region, with d = 0.2, 
the bulk of the work is done by SSNR and little 
by SNR. A value of S in the range 0.15 to 0.2 is 
indicated, with 0.15 being more physical. In the 
Outer Galaxy , the likely range is 0.0 to 0.15. 

Certainly, the values of / shown in Figure 4 
have big systematic and statistical errors due 
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mainly to the simplified model of clustering and 
big fluctuations caused by the stochastic nature 
of SN explosions. The estimates of these errors 
are indicated in Figure 4. 

5. CONCLUSIONS 

The fluctuations in the intensity of gamma rays 
above 1 GeV in and near the Galactic Plane give 
information about the mode of production of cos- 
mic rays in their sources, and their mode of prop- 
agation. 

We have adopted alternative models of both 
sources ( supernovae of a variety of energies and 
SN of unique energy but with varying degrees 
of 'clustering' ) and diffusion ( our 'standard' 
anomalous diffusion, and normal ). For 'nor- 
mal diffusion', although at D ~ 5.5 and 11.5 kpc 
it is possible to achieve 'reasonable' results, for 
7 < D < 11 kpc ( the region nearest the sun ) 
this is not possible. Whatever the value of 5, 
the fractions are unreasonably high. Perhaps this 
argument can be used to point against 'normal' 
diffusion being applicable for CR of the energies 
in question ? 

Bearing in mind our general arguments that 
clustering is more common in the Inner Galaxy it 
would appear that S ~ 0.15 is favoured ( Figure 
4b ). We would then have a constant fraction of 
SN of the SSN variety, with none of them clus- 
tered, but the normal SN clustered only in the 
solar vicinity and the Inner Galaxy. 

It seems that the solution put forward is astro- 
physically reasonable, for the following reasons. 
Concerning standard SNR, these come from stars 
of modest mass produced in Molecular Clouds. 
These clouds are composed of individual clumps 
in which star production ( and subsequent SN ) 
occurs. In the Inner Galaxy, where the overall 
MC masses are bigger [T2J there will be more clus- 
ters than in the Outer. For SSN, the fractional 
number per clump of the required massive pro- 
genitor stars is probably the same in the Inner 
and Outer Galaxy. In view of continuing prob- 
lems with theories of star formation [ T.J.Millar, 
private communication ] even this conclusion can- 
not be regarded as completely firm. 

As for the UHECR connection Figure 5 shows 



the energy spectra of CR created by SSNR. The 
distribution of their explosion energies and max- 
imum energies of accelerated CR, taken from [3] 
and used in our analysis, gives the possibility to 
accelerate particles up to ultra-high energies of 
about 10 9 GeV. So our preferred explanation of 
large fluctuations for the gamma-ray intensity has 
an additional advantage to include the possible 
origin of UHECR. 
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Figure 5 . Energy spectra of CR accelerated by two types 
of SN - SNIbc and SNIIn: (a) the summed spectra of 10000 
SNIbc and 5000 SNIIn; (b) the mean spectrum and its ir- 
regularity for samples shown in (a) ; The irregularity is the 
standard deviation from the mean CR intensity for dif- 
ferent samples of the spectrum indicated by vertical lines. 
The curve obs(P) is our estimate of the proton spectrum 
( or the rigidity spectrum for heavier nuclei ) needed to fit 
the experimental data \ '.'> . 
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Table 1 





A£ 


No. of entries 


Mean 


a 


Inner 


1° 


161 


0.00295 


0.087 




2° 


80 


0.00225 


0.086 




4° 


40 


0.00175 


0.081 




8° 


20 


0.000500 


0.074 


Outer 


1° 


180 


-0.00272 


0.098 




2° 


90 


-0.00156 


0.096 




4° 


45 


-0.00167 


0.091 




8° 


23 


+0.00022 


0.081 



Table 1. Results for the distribution of the differences 
between the EGRET intensities and the 9- degree polyno- 
mial fit: nos. of entries, means and RMS (a-) values. All 
values relate to the logarithm of the intensity. Latitude 
range, \b\ < 2.5°, E 1 > 1 GeV. The histograms are given 
in Figure 2. 
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